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I Tevatron data, as yet there is no unambiguous evidence for other than standard next-to-leading- 

' order GLAP evolution. The framework of dipole models and saturation of parton densities is 

o ■ 
o 

I of the data, so do other, more conventional, explanations. 

5^ ■ 1 Introduction 



Abstract 

Low-x physics is reviewed, with particular emphasis on searches for deviations from GLAP evolu- 
tion of the parton densities. Although there are several intriguing indications, both in HERA and 



examined and confronted with the data. Although such models give a good qualitative description 



The scattering of energetic 'simple' particles from an unknown target to elucidate its 



H ' structure is an experimental technique with a long and distinguished history. Such seat- 

ed ' 

tering experiments have revolutionised our view of the microscopic world; the prototype, 
and most famous, is the scattering of alpha particles from a thin gold foil carried out by 
a previous president of this society. Lord Rutherford, working with Geiger and Marsden 
in Manchester in 1909. This experiment led to the concept of the nuclear atom 
A similar revolution occurred after the SLAG experiments, as discussed by Taylor at 
this meeting These led to the general acceptance of the concept of quarks as actual 
constituents of the nucleon, rather than as mathematical abstractions. 

With the advent of the HERA electron-proton collider, the explorable phase space in the 
kinematic invariants (the virtuality of the exchanged virtual photon) and x (the frac- 
tional momentum of the parton involved in the scattering) has increased by approximately 
three orders of magnitude in each variable (see figure |l]). This increase in kinematic range 
has opened up a new branch of studies, which can be described generically as 'low-x 
physics'. As will be seen in this talk, the study of this kinematic region, which for conve- 
nience will be defined hj x < 10~^, is full of interest and has already led to many advances 
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Figure 1: The kinematic plane in x and for all experiments, probing the 
parton distribution of the proton. The region of interest in this talk is indicated to 
the left of the vertical line. The other line at close to 45° indicates the approximate 
kinematic limit at HERA 

in tlie understanding of tlie theory of tlie strong interaction, Quantum Chromo dynamics 
(QCD). 

Although the study of diffractive processes is intimately linked to many aspects of low-x 
physics, constaints of time mean that it is not covered in this talk. The contribution by 
J. Dainton |Q to this meeting touches upon diffraction to some degree. 

1.1 Theoretical background 

The interest in low-x physics is that particles with small x are the result of a large number 
of QCD branching processes. The behaviour of partons at low x thus reflects the dynamics 
of QCD and allows the behaviour of its couplings and interactions to be probed over a 
large range in the kinematic variables. In particular, the evolution of the number of 
partons as a function of x and will be sensitive, depending on the kinematic range, to 
the various approximations that describe QCD evolution. 

Several of the contributions to this meeting have discussed the subject of QCD evolution 
in some depth ||5|,^, and there are of course many excellent overviews of the subject 0, so 
that is appropriate here only to give a very brief summary of the most important points 
relevant at low x. 

One of the most important properties of QCD, without which its usefulness as a theory 
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would be extremely limited, is that of factorisation. This states that hard processes can 
be regarded as a convolution of a 'sub-process' cross section that can be calculated in 
terms of point-like interactions together with the probability to find the participating 
particles in the target and in the probe. The subsequent hadronisation of the participants 
in the hard collision, together with the target and probe remnants, can be regarded as an 
approximately independent process. Thus the cross section can be written schematically 
as: 

(y^ft®U®a (1) 

where a is the sub-process cross section and /t, fp are the parton distribution functions for 
the target and probe, respectively. One of the most important results of the factorisation 
hypothesis is that the parton distribution functions (PDFs) measured in one process can 
be used in the cross-section determination for a completely different process. Furthermore, 
QCD provides the tools by which to extrapolate from the PDFs measured at one scale to 
very different scales. 

Specialising now to deep inelastic scattering (DIS), the PDF for the highly virtual photon 
can normally be considered to be a 5 function, so that equation |l] becomes: 

a ~ / (g) a (2) 

where / now represents the PDF of the proton. It is conventional to assume that / 
satisfies the schematic equation: 

- "^^^'^ {f^V) (3) 



where /i represents the renormalisation scale and P is a 'splitting function' that describes 
the probability of a given parton splitting into two others. This equation is known as 
the (Dokshitzer)-Gribov-Lipatov-Altarelli-Parisi equation IsHTll. There are four distinct 
Altarelli-Parisis (AP) splitting functions representing the 4 possible 1-^2 splittings 
and referred to as Pgg, Pgq, Pqg and Pgg. The calculation of the splitting functions in 
perturbative QCD in equation ^ requires approximations, both in order of terms which 
can be taken into account as well as the most important kinematic variables. The generic 
form for the splitting functions can be shown to be : 



= E -4!:' {'■>(;)} 

n=0 m=0 \ / ) 



m 



(4) 



where as is the strong coupling constant, V (x) are the x-finite parts of the AP splitting 
functions and Am^ are numerical coefficients that can be calculated, at least in principle. 



3 



for each splitting function. The AP sphtting functions sum over terms proportional to 
(aglnQ^)" in the perturbative expansion. Thus, for example, the term in equation ^ 
with n = m = when added to V^^^ (x) corresponds to leading order in so-called Gribov- 
Lipatov-Altarelli-Parisi (GLAP) evolution. In some kinematic regions, and in particular 
at low X, it must become essential to sum leading terms in In 1/x independent of the value 
of In Q^. These terms in some sense correspond to corrections taking into account so-called 
Balitsky-Fadin-Kuraev-Lipatov |T^-|T5[ (BFKL) evolution, which governs the evolution in 
X at fixed Q^. As x falls, this must at some point drive parton evolution. One of the 
continuing themes of low-x physics, as will be discussed in sections ^ and ^, is the search 
for experimental effects that can be unambiguously attributed to BFKL evolution. 

Figure shows the Inl/x - InQ^ plane at HERA, together with schematic indications of 
the directions in which GLAP and BFKL dominate the evolution of parton distributions. 
The third direction on the figure, labelled 'CCFM', refers to an approach to an integrated 
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Figure 2: Schematic diagram showing different regions of the In l/x and InQ^ 
plane and the evolution equations expected to hold therein. The line marked 'sat- 
uration' represents the boundary between GLAP evolution and evolution governed 
by the GLR equation. The 'size' of partons is also indicated in differing kinematic 
regions. 

evolution containing both the leading GLAP and BFKL terms to equal order developed 
by Ciafaloni, Catani, Fiorani and Marchesini [p!6HT8[|. Also indicated on the figure are 
schematic indications of both the 'size' and density of partons in the proton in different 
kinematic regions. The transverse size of the partons which can be resolved by a probe 
with virtuality is proportional to 1/Q, so that the area of the partonic 'dots' in figure]^ 
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falls as rises. For particular combinations of parton size and density, the proton will 
eventually become 'black' to probes, or, equivalently, the component gluons will become 
so dense that they will begin to recombine. The dotted line labelled 'Critical line - GLR' 
refers to the boundary beyond which it is expected that such parton saturation effects 
will become important, i.e. the region in which partons become so densely crowded 
that interactions between them reduce the growth in parton density predicted by the 
linear GLAP and BFKL evolution equations. The parton evolution in this region can be 
described by the Gribov-Levin-Riskin [19,^] equation, which explicitly takes into account 



an absorptive term in the gluon evolution equation. Naively, it can be assumed ||2T| that 
the gluons inside the proton each occupies on average a transverse area of 7rQ~^ so that the 
total transverse area occupied by gluons is proportional to the number density multiplied 
by this area, i.e. '7iQ~'^xg{x,Q'^). Since, as will be discussed later, the gluon density 
increases quickly as x falls, and the gluon 'size' increases as Q~^, in the region in which 
both X and are small, saturation effects ought to become important. This should occur 
when the size occupied by the partons becomes similar to the size of the proton: 

xg{x,Q')^ = 7rR' (5) 
where R is the radius of the proton, (~ 1 fm ~ 5 GeV~^). The measured values of 



xg{x,Q ) imply that saturation ought to be observable at HERA |22| at low x and Q , 



although the values of which satisfy equation |^ are sufficiently small that possible non- 
perturbative and higher-twist effects certainly complicate the situation. Of course, it is 
also possible that the assumption of homogenous gluon density is incorrect; for example, 
the gluon density may be larger in the close vicinity of the valence quarks, giving rise 
to so-called 'hot spots' [^], which could lead to saturation being observable at smaller 
distances and thereby larger Q^. The concepts of 'shadowing' or saturation have been 
discussed now for many years [p0| , |2^ ,|2^- p^j . As will be seen in section HERA does 



indeed provide data of relevance to such discussions. 



2 The Structure Function Data 

In this section, the most recent structure function data from ZEUS and HI are presented 
and discussed. After some initial definitions of kinematic variables and the structure 
functions relevant at low x, the data on F2 are shown and indirect methods of extracting 
the longitudinal structure function, Fl are discussed. 
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2.1 Kinematics and structure function formulae 



The scattering of a lepton from a proton at sufficiently large can be viewed as the elastic 
scattering of the lepton from a quark or antiquark inside the proton. As such the process 
can be fully described by two relativistic invariants. If the initial (final) four-momentum 
of the lepton is k{k'), the initial four-momentum of the proton is P, the fraction of the 
proton's momentum carried by the struck quark is x and the final four-momentum of the 
hadronic system is P', the following invariants may be constructed: 

s = {P + kf (6) 
= -q' = -{k'-kf (7) 

= {P'f = {P + qf (9) 
Energy- momentum conservation implies that: 

so that, ignoring the masses of the lepton and proton: 

y = - (11) 

sx 

= Q^l—^^^ (12) 

X X 

where the approximate relationship will in general be sufficiently accurate for the values 
of X of interest in this talk. Since DIS at a given s can be specified by any two of these 
invariants, the most convenient may be chosen, normally x and Q^. 



Equation shows the general form for the spin-averaged neutral current differential cross 
section in terms of the structure functions Fi, F2 and F^: 

(fa 2tiq^ 



[2xy''Fi + 2(1 - y)F2 



dxdQ^ xQ^ 

± {1 - {1 - yfjxFs] (13) 

where the -|- sign in the ± term is taken for e~ and the — sign for interactions. The 
structure functions are products of quark distribution functions and the couplings of the 
current mediating the interaction. They are in general functions of the two invariants 
required to describe the interaction. 

To leading order in the QCD-improved parton model, in which quarks are massless, have 
spin 1 and in which they develop no pr, the Callan-Gross relation |[53|] : 

2xFi{x) = F2{x) (14) 
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is satisfied. At the next order, px must be taken into account and this relation is violated. 
This is usually quantified by defining a longitudinal structure function, Fl, such that 
Fl = F2 — 2xFi. Substituting into equation [1^ gives: 



dxdQ'^ 



[Y+-F2ix,Q' 



xQ* 

y'FL{x,Q^) + Y.-xFs{x,Q^)] 



(15) 



where Y± are kinematic factors given by: 

F± = 1 ± (1 



(16) 

At low X, in general <^ Af|, xF^ vanishes and equation reduces to that for photon 
exchange. In the rest of this talk, electroweak effects will be neglected. 

In general, the form of the structure functions beyond leading order depends on the renor- 
malisation and factorisation scheme used. In the so-called 'DIS' scheme, the logarithmic 
singularity produced by collinear gluon emission is absorbed into the definition of the 
quark distribution, so that the structure functions have a particularly simple form and 
can be expressed to all orders as 

F2{x,Q') = 

i=u,d,s,c,b 



J2 A{Q')[xqi{x,Q')+xq,{x,Q')] 



(17) 



The parton distributions qi{x, Q"^) and gj(x, Q^) refer to quarks and antiquarks of type i. 
The quantities Ai{Q'^) are given by the square of the electric charge of quark or antiquark 



i. However, in the MS scheme, the form of F2 changes with order in QCD. In LO QCD 
it has the form: 

1 



dy f X 



y \y 



5{l- 



y 



{yq^{y. Q') + yq.{y. Q')) + <| ^Cf [-]) yg{y, Q 



27r 



X 



(18) 



where g{x,Q'^) is the gluon density in the proton, ^^(Q^) is the QCD running coupling 
constant and Cq{x) and Cg{x) are scheme-dependent 'coefficient functions'. 

In contrast, the longitudinal structure function contains no collinear divergence at first- 
order in QCD so that: 

1 



Fl{x.Q' 



+ 



2 /'^^^ 



i=u,d,s,c,b 
2 



y \y 



X 



(19) 
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independent of the factorisation scheme employed. 

Taking account of radiative corrections via the term 5^, equation [TB becomes: 

^ = %-i^+5r)-K-F,M^)-y^F,{.,Q^)] (20) 

A useful quantity known as the 'reduced cross section' can be defined from equation [2^ 
by taking the kinematic factors to the left-hand side, i.e.: 

xQ^ (Pa 



27ra2y+ ■ (1 + 5) dxdQ"^ 



= F2{x,Q')-^Fl{x,Q') (21) 
+ 

and, provided y is small, ar is to a good approximation equal to F2. 

An alternative formalism to describe DIS interactions at low x in terms of total virtual 
photon-proton cross sections is particularly useful when discussing the low-Q^ region and 
the transition to real photoproduction. The total cross section can be written as the sum 
of the cross sections for transversely and longitudinally polarised virtual photons: 

atr{W\Q')^aT + aL (22) 



where x and W are related using equation 12. This leads directly to expressions for F2 



and Fl in terms of virtual photon-proton cross sections: 

atriW',Q')r^ F2 =^{aT + aL) (23) 
Fl =j^^l (24) 

2.2 The F2 data at 'medium' 

In this section the most recent F2 data at 'medium' values of from the two HERA 
experiments are discussed. The term 'medium Q^' used here is essentially a definition 
related to the characteristics of the HI and ZEUS detectors; the term covers the struc- 
ture function measurements in which the minimum scattered electron angle (and hence 
the minimum Q2) that can be measured is determined by the dimensions of the main 
calorimeters of the two experiments. This is in contrast to the 'low-Q^' region, where, as 
discussed in section |2.3| , the minimum range is determined by the geometric acceptance 
of small, special-purpose detectors placed downstream in the electron-beam direction. 

The measurement of F2 is a very complex and painstaking effort and has been often 



described before [p^p5|. Since both the HI and ZEUS detectors are sufficiently hermetic. 
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the two invariants required to define the event kinematics fully can be reconstructed from 
measurements on the electron, on the hadronic final state corresponding to the struck 
quark, or on mixtures of the two. The optimal method depends on the kinematic region 
of interest and on the properties and resolutions of the detectors. The size of the radiative 
corrections is also very dependent on the reconstruction method employed. The basic 
measurements made are the angles and energies of the electron and current jet. Many 
methods have been used by the two experiments, of which probably the most important 
for the measurement of F2 are: the Electron Method, in which the energy and the angle 



of the scattered electron are used; the Double Angle method [0], in which the angles of 
the scattered electron and the current jet are used; the S method which uses the 
current jet and electron energies and the electron angle, and the Pt method |58], which 
uses the Double Angle method with the additional constraint of pt balance. 

In order to determine F2 the following steps are carried out. First, a sample of DIS 
events is selected, basically by requiring an identified electron in the detector. The kine- 
matic variables are reconstructed using one of the methods discussed above. The data are 
binned in x and with bin sizes determined by detector resolution, statistics, migration 
in and out of the bin due to the finite resolution of the experiments, etc. Estimates of 
the background in each bin are made and statistically subtracted. The background in 
the low-x region is dominated by photoproduction processes in which a fake electron is 
reconstructed because of confusion with hadronic debris from photoproduction interac- 
tions, which of course have a much higher cross section than DIS processes. The data 
are corrected for acceptance, radiative effects and migration via Monte Carlo simulation. 
Multiplying the corrected number of events by the appropriate kinematic factors shown 
in equation ^ and subtracting an estimate for gives values of F2 and hence the quark 
and antiquark densities in the chosen bins. There are great gains in physics terms to 
be made by pushing the precision of these measurements to the limits. Since for much 
of the available phase space systematic effects are dominant, this requires progressively 
better understanding of the detectors involved and simulation of their response to parts 
per mille. 

Figure |^ shows the preliminary HI measurement of the reduced cross section (see equa- 
tion ^) in bins of as a function of x. Also shown are data from the fixed-target 
experiments NMC and BCDMS |60|,|T|. The bins from Q"^ = 1.5 GeV^ to 150 GeV^ 



are shown. In the relatively small region of overlap, there is good agreement between the 
HI and fixed-target data. The most obvious characteristic of the data is the steep rise 
of F2 at low X. The curve shown on the figure is the result of an NLO QCD fit by the 
HI collaboration to this data together with earlier measurements at higher x from NMC. 
The NLO QCD fit [^, based on CLAP evolution, uses three light flavours with charm 
added via the boson-gluon fusion process and uses as(M|) = 0.118. It gives an excellent 
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Figure 3: The preliminary HI data on the reduced cross section (proportional 
to F2) from the 1996-97 data-taking period. Also shown are points from the fixed- 
target experiments NMC (triangles) and BDCMS (squares). The solid curve shows 
the NLO QCD fit carried out by HI, while the dotted curve visible at the lowest x 
corresponds to the expectation for Fl = 0, as discussed in the text. 



fit to the data over the full kinematic range. The quality of the QCD fit permits the 
conclusion that the rise of F2 at low x is unambiguously associated with a dramatic rise 
in the gluon distribution. Also shown in figure |] as the dashed line is the expectation of 
the fit for F2 alone. As can be seen, there is a small departure from the measured value of 
(Tr, implying that in this kinematic region, the effect of Fl begins to become perceptible; 
this is discussed further in section E^. 
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Figure § shows the prehminary 1996-97 ZEUS data, plotted as the F2 structure function 



ZEUS Preliminary 96/97 
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Figure 4: The lowest bins of the preliminary ZEUS data on F2 from the 
1996-97 data-taking period. The preliminary data are in good agreement with the 
published data (shown as open circles). Also shown are points from fixed-target 
experiments. The solid curve shows the ZEUS NLO QCD fit, while the dashed line 
shows the CTEQ4D curve and the dotted line that from MRST99. 



as a function of x in bins with the QCD prediction subtracted. Only the lowest 
bins are shown; F2 has been determined up to = 30,000 GeV^. Fixed target data from 
BCDMS [HP, E665 |63|, NMC pf, and SLAG H are also shown. The data agree 
well with the HI data and show the same dominant feature of a very steep rise at low x. 
The data are also very well described by the NLO QCD global fit to parton distributions 
of CTEQ4D m and MRST99 ||. 

Figure ^ shows the HI data together with data from NMC and BCDMS, now plotted in 
X bins as a function of Q^. The data cover approximately five orders of magnitude in 
both X and Q^. At high x, approximate scaling in can be clearly observed. As x falls, 
deviations from scaling become stronger and stronger. The lines on the figure are the 
result of the NLO QCD fit, which can again be seen to give an excellent description of 
the data. 



The HI data for F2 at low x in x bins as a function of In is shown in figure |^. It can 
clearly be seen that the data are not linear in InQ^. In fact they fit well to a second-order 
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Figure 5: Preliminary HI data on F2 from 1996-97 and published data from 
1994-97 in bins of x as a function of Q"^. Also plotted are fixed-target data from 
SLAC, NMC and BCDMS. Each x bin is offset by the amount indicated in the 
legend for ease of visibility. The curves show the HI NLO QCD fit to the data. 
The flat 'scaling ' regime at high x gives way at lower x to steep scale breaking due 
to gluon radiation. 



polynomial of the form 

F2 = A{x) + B{x) lnQ2 + C(x)(lnQ2)2 (25) 



and this polynomial fit is almost indistinguishable from the HI NLO QCD fit, except at 
the highest Q^. 

The most convenient and useful way to parameterise the deviations of the data from 
scaling is to examine the logarithmic derivative, dF2/d\nQ^ , which in leading-order QCD 
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Figure 6: Preliminary HI data on F2 from 1996-97 in bins of x as a function 
of . The solid line shows a fit to the form of equation ^ while the dashed line 
shows the HI NLO QCD fit. 



is directly proportional to the gluon density at twice the x of the derivative^ 

dF2 la. 



(91ng2 



97r 



■xg{x, 



(26) 



Having fit the data to the form of equation it is straight-forward to obtain the log- 
arithmic derivative. Since this is proportional to the LO gluon density, it is clear that 
a precision measurement of the scaling violations can be used directly to determine the 
gluon distribution. Figure |^ shows such determinations from both HI and ZEUS. The HI 
determination comes from the NLO QCD fit referred to above whereas that from ZEUS 
comes from an NLO fit to published data The steep rise in the gluon density as x 
falls is apparent. Also noticeable, particularly in the ZEUS determination, is that this rise 
becomes weaker and weaker as falls. Indeed, for = 1 GeV^, the gluon density falls 
below that of the singlet quark structure function and is essentially compatible with zero. 
This seems to contradict the 'standard' picture in which the rise of F2, which is of course 
only directly sensitive to the density of charged partons, is driven by the quark-antiquark 
pairs produced from gluons. However, these interesting effects only become obvious at 



^ It should be noted that, although a good approximation at LO, this relation becomes increasingly less 
valid at higher orders. Although a very useful qualitative relationship, it should therefore be used with 
circumspection. 
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Figure 7: The left-hand plot shows preliminary HI data on determination of the 
gluon from F2 data from 1996-97. The gluon density is shown for three values of 
Q^; the central band shows the statistical uncertainty while the outer bands show 
the systematic and the theoretical uncertainties added in quadrature. The right- 
hand side plot shows ZEUS data on determination of the gluon from F2 data from 
1995. The gluon density is shown for three values of ; all data down to > 1 
GeV^ is included in the NLO QCD fit. The bands show all uncertainties added in 
quadrature. The light shaded band is the gluon density while the dark shaded band 
shows the sea. 



very low values of Q^, and it is not clear that it makes sense to talk about 'gluon densities' 
at these low valuesQ Nevertheless, it is certainly true that the NLO QCD fits them- 
selves seem to give a perfectly satisfactory description of the general features of the data 
down to these low values of Q^. Much more will be said on this subject in section O. 



In addition to the high-precision data on the fully inclusive F2, the ZEUS collaboration 
has also presented data on semi-inclusive DIS in which a charm quark or antiquark is 
involved in the hard scatter [^]. Figure ^ shows the F^^^ data in x bins as a function of 
Q^. A qualitatively similar pattern of scaling violations to that in the fully inclusive F2 
can be seen; however, the scaling violations seem, within the relatively large errors, to be 
stronger than in the inclusive case and to set in rather earlier. While part of this effect 
can be attributed to the effect of the charm-quark mass, it is also to be expected since 
the dominant process in DIS charm production is boson-gluon fusion, which is entirely 
driven by the gluon density in the proton. Once again, it can be seen that the NLO QCD 
fit gives an excellent description of the data. Figure ^ shows the ratio of the charm over 
the inclusive F2. At small x the ratio flattens, implying that the charm and inclusive 
structure functions grow at the same rate, as is to be expected if both are dominated by 



^ It was interesting to note Dokshitzer's comment in the discussion sessions at this meeting that in fact 
it docs make sense to discuss gluon distributions at such low values of Q^. 
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Figure 8: ZEUS data on the charm structure function, in bins of x as a function 
ofQ^. The band shows the ZEUS NLO QCD fit. 

the gluon in this region. For low values of x, the ratio falls at fixed x as falls. This 
is consistent with the observation discussed above that as falls the gluon density at 
fixed X also falls. 



The ZEUS and HI detectors are not perfectly hermetic, since it is clearly necessary to 
allow the beams to enter and leave the apparatus. Thus the 'beam-hole' limits the angular 
acceptance of the detectors both at very forward and very backward directions. In the 
very backward direction (small lepton-scattering angles) this limits the values that 
can be accurately measured to around ~ 2 GeV^. In order to access smaller (and 
thereby smaller x) , the geometrical acceptance of the detectors must be extended in some 
way. There are two main ways in which this has been achieved. The first is to shift 
the interaction vertex in the direction of the proton beam, typically by of order 60 cm, 
so that the electron has further to travel before it strikes the rear calorimeters. This 
means that the geometrical edge of the detector now corresponds to a smaller scattering 
angle, and hence lower can be accepted. The other method is to install small, high 
precision, detectors further upstream in the electron beam direction which can thereby 
detect much smaller scattering angles than the main detectors. Both ZEUS and HI have 
such detectors, although so far only ZEUS have published results. 

ZEUS published some time ago results using their Beam Pipe Calorimeter (BPC) ||69|| , 



2.3 The F2 data at 'low' 
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ZEUS 1996-97 




Figure 9: ZEUS data on the charm structure function, plotted as a ratio of the 
inclusive structure function, F2 in bins of as a function of x. The band shows 
the ZEUS NLO Q CD fit. 

which is a small tungsten-scintillator sampling calorimeter placed 2.94 m away from the 
interaction point in the electron beam direction. In 1997, two silicon-microstrip detector 
planes were added in front of the calorimeter in order to improve the position resolution. 
ZEUS has recently published the final results from this BPC/BPT combination |7^, which 



extend the measurement of F2 down to a; ~ 6 ■ 10~^ and ~ 0.045 GeV^. Figure 10 
shows the final ZEUS data in bins of as a function of x. The new data match well 
with the previous ZEUS BPC data, as well as with that from other experiments in the 
overlap region. However, the extrapolation of the ZEUS Regge fit (see below) into the 
fixed target regime is generally of order 15% above this data. 

The solid curve labelled 'ZEUS Regge fit' on figure [T^ shows the result of a fit to the form: 

where A-^, A-p and Mq are constants and ajR and ap are the Reggeon and Pomeron 
intercepts, respectively. This phenomenological parameterisation is based on the combi- 



nation of a simplified version of the generalised vector meson dominance model for 
the description of the dependence and Regge theory ||7^ for the description of the x 
dependence of F2. Regge Theory is most applicable to the description of cross sections at 
asypmtotic energy. Equations ^ and ^ which relate F2 to cross sections evaluated at 
energies proportional to imply that Regge theory should describe the very low-x data 
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Figure 10: ZEUS BPT data on F2 in bins of as a function of x. Also shown 
are earlier ZEUS data as well as data from HI and E665. The solid line shows 
the results of the 'ZEUS Regge fit' to the form of equation while the dotted line 
shows the result of the ZEUS NLO QCD fit. 



well. This is borne out by figure |10], where the ZEUS Regge fit gives a good description of 
the data up to ~ 1 GeV^. Above this Q^, however, the Regge description rapidly fails, 
whereas the ZEUS NLO QCD fit, shown for > 1 GeV^, is an excellent description of 
the data from here to the highest . 



Figure |TT| shows the ZEUS F2 data in bins of constant y as a function of InQ^. For 
Q2>1 GeV^, the data are roughly independent of Q^, whereas at lower they fall 
rapidly, approaching the fall-off that would be expected in the limit — > from 
conservation of the electromagnetic current. Whether this dependence indicates that 
this limit has already been reached, or whether other effects, for instance saturation, are 



responsible, will be discussed in section 4.4 
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Figure 11: ZEUS BPT data on F2 in bins of y as a function ofQ"^. Also shown 
are earlier ZEUS data as well as data from HI and E665. The solid line shows 



the results of the 'ZEUS Regge fit' to the form of equation while the dotted line 
shows the result of the ZEUS NLO QCD fit. 



2.4 The Fr structure function 



As discussed in section the differential cross section for DIS at low x depends on two 
structure functions, F2 and F^. Since in principle both F2 and Fl are unknown functions 
that depend on x and Q^, the only way in which they can be separately determined is 
to measure the differential cross section at fixed x, and at different values of y, since 
as shown in equation |15], the effect of F^ is weighted by y^ whereas F2 is weighted by 
1 + (1 — y)"^. However, since = sxy, fixed x and implies taking measurements 
at different values of s. This can certainly in principle be accomplished by reducing the 
beam energies in HERA. However, the practical difficulties for the experiments and the 
accelerator inherent in reducing either the proton or electron beam energy, or both, by a 
factor sufficient to permit an accurate measurement of Fl mean that it has not to date 
been attempted. An alternative way to achieve the same end is to isolate those events 
in which the incoming lepton radiates a hard photon in advance of the deep inelastic 
scattering, thereby reducing the effective collision energy. Unfortunately, the acceptance 
of the luminosity taggers typically used to detect such photons is sufficiently small and 
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understanding the acceptance sufficiently difficult that, although both experiments are 
working on the analysis, neither has as yet produced results. 

In the absence of any direct determination, the HI experiment has utilised its ability to 
detect events at very large values of y in order to carry out an indirect measurement of 
Fl. The determinations of F2 rely on the fact that most of the measurements are made 
at values of y sufficiently small that the effects of Fl are negligible; those at higher y 
usually have an estimate of F^, which according to QCD is in any case normally a small 
fraction of F2, subtracted off. The HI collaboration inverts this procedure by isolating 
kinematic regions in which the contribution of F^ is maximised and then subtracts off the 
QCD prediction of F2 measured at lower y. 

As remarked earlier, figure ^ shows the reduced cross section, defined by equation ^ in 
which the contribution of Fl can be seen at the lowest x (which, for fixed Q^, corresponds 
to the highest y) as the difference between the full QCD fit and that with Fl set to zero. 
Thus, Fl can be estimated from the following relationship: 

W - ( pQCDfit XQ"^ „ \ ^+ 

- V' ~ 2^^^J ■ ^ ^^^^ 



An alternative method used by HI employs the derivatives of the reduced cross section 
with respect to Iny, thereby making rather different QCD assumptions. Differentiating 
equation ^ leads to the following expression: 

dar dF2 2y\2-y) OFl 



dlny d\ny Yf dlny 

which leads to greater sensitivity to Fl via the stronger y dependence at the cost of 
involving derivatives of (Jr,F2 and Fl, the quantity to be measured. It is instructive to 
consider various restrictions: 

• Small y - here dar/dlny ~ dF2/dlny. For low x, F2 can be well approximated by: 

F2 oc x"^ oc y^ (30) 



so that: 



which can be expanded as: 

dF2 



^""^ xy' (31) 



dlny 



oc 



Xe^^'^y ^ X{l + X\ny...) (32) 



d\ny 

provided Xlny is small. From this it is clear that dar/dlny is linear in Iny; 
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Figure 12: Preliminary HI data on the derivative of the reduced cross section 
with respect to \ny in bins. The solid curve shows the result of an NLO QCD 
fit with the value of Fl as predicted by QCD from the measured F2. The dotted line 
shows the same fit with Fl = and the dashed line shows the fit with Fl = F2. 

• Fl = - for all y, dar/dlny is linear in \ny for the same reason as above; 

• Fl ^ and large y - dcr/dlny is non-linear in Iny and the deviations are proportional 
to Fl and its logarithmic derivative; 

• large at small y - this implies x becoming larger so that at some point the ap- 
proximation of equation ^ starts to fail and therefore there are deviations from non- 
linearity. 

All of these features can be seen in the preliminary HI data of figure At the largest 
values of y, the deviation from linearity implies that Fl is non-zero. Although it is 
in principle possible to solve the differential equation for Fl implied by equation ^ in 
practice the data are insufficiently precise, so that the value of the derivative is taken from 
the QCD fit. Variations in this are included in the systematic error. Also in principle it 
is possible to iterate the Fl estimated in this way with that assumed in the measurement 
of F2; once again, the precision of the data does not permit this and in any case the 
correlation would be very large. 

The HI collaboration have used the first method discussed above to estimate Fl for 
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Figure 13: Preliminary HI estimate of F^. The Fl values obtained are plotted in 
bins as a function of x. Also shown are earlier bins at higher x from the SLAC 
and NMC experiments. The solid line is the prediction of the HI NLO QCD fit. 



Q > 10 GeV and the second for smaller Q . The results are shown in figure 13, together 



with earlier determinations from SLAC ||6J], NMC |^ and BCDMS [^. The curve is 
the result of an NLO QCD fit to the HI data deriving from the F2 determination, i.e. by 
deriving the gluon and quark distributions from scaling violations and then calculating 
Fl using a QCD formula such as equation |l^. The QCD prediction is in good agreement 
with the HI estimate. 

In summary, although the indirect determinations of F^ are both interesting and impor- 
tant, there is no substitute for a direct measurement. Since after the HERA upgrade the 
lowest regime will no longer be accessible because of the new final-focus quadrupoles 
that close off the small scattering angle aperture, such data will presumably have to come 
from hard initial-state radiation events. First results from HI and ZEUS are eagerly 
awaited. 



3 Other probes of QCD dynamics at small x 

Despite the very high precision and very large kinematic range of the structure-function 
data shown in the previous section, there was no obvious sign of any deviation from NLO 
CLAP evolution (although see section [4.4|) . Indeed, there are very good reasons why this 



should be so, even if BFKL dynamics were important in this kinematic range. It is 
generally agreed that the chances of observing any deviation from CLAP evolution are 
greatly enhanced by examining certain exclusive processes in particular corners of phase 
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space. Several such processes are examined in this section: the production of forward jets 
and forward vr" at HERA and the study of dijets with a large separation in rapidity at 
the Tevatron. 



3.1 Forward jet production at HERA 

One of the most marked characteristics of GLAP evolution is that successive parton 
branchings are strongly ordered in k^, the square of the transverse momentum of the 
parton. In the case of BFKL, there is essentially a random walk in k^, while strong 
ordering occurs in 1/x. These observations immediately imply the corner of phase space 
that is most likely to exhibit BFKL effects: small x, since this will enhance the importance 
of the ln(l/x) terms in the pQCD expansion, and pt of the struck quark ~ Q'^, which will 
strongly suppress GLAP evolution because of the strong kr ordering of successive parton 
branchings between the virtual photon and the struck quark |75[. The kinematic 



properties of the struck quark can be reconstructed in several ways, of which the most 
usual is to tag a high-energy jet. The kinematic requirements discussed above imply 
selecting events with low x and therefore low Q^, so that the balancing jet with comparable 
Pt will be in the very forward direction. 

The HI collaboration presented preliminary results on forward jet production at the 
DIS2000 conference in Liverpool |j76|, in which they isolated jets using the inclusive k^ jet 



algorithm in two pseudorapidity regions: 'central', defined as 0.5 < t] < 1.5, and 'forward', 
defined as 1.5 < ?7 < 2.8. Jets were selected with fractional energy x\et = p '"^ > 0.035. 

-^pbcam 

The differential cross-section da-^ct/dx for the 'forward' and 'central' regions is shown in 
figure |14[ Whereas NLO QCD gives a reasonable description of the data in the central 
region, it clearly falls below the data in the forward region at the lowest values of x. 
Although the hadronisation corrections are largest in the forward direction, they are 
insufficient to explain the discrepancy. The variation in the NLO QCD prediction by 
varying the scale by a factor of two in each direction is also large, but again insufficient 
to explain the shortfall. However, is used as the hard scale in these calculations, and 
it is entirely unclear whether this, or Q^, is the appropriate scale. Figure ^ shows a 
comparison of the calculations assuming and to give the hard scale in the NLO 
calculation. It can be seen that not only is the discrepancy reduced when is used, the 
uncertainty in varying the scale from Q"^ /2 to 2Q^ is sufficient to give agreement with the 
data. Given this theoretical uncertainty, it is clear that no firm conclusion can be drawn 
on the presence of non-GLAP evolution in forward jet production. 
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Figure 14: Preliminary HI data on production of inclusive jets, found with the 
inclusive kT jet algorithm, in two bins of rj as a function of x. The dotted line 
shows the result of an LO QCD fit, whereas the dotted line within the shaded hand 
shows the NLO QCD fit at a scale of ^\ = with an estimate of the theoretical 
uncertainty. The central plots show the effect of hadronisation corrections, while 
the difference of the data from the NLO predictions is shown at the bottom. 



3.2 Forward tt" production at HERA 



Since the size of any BFKL effect will certainly strongly increase as x falls, processes 
that can reach lower values of x than possible for forward jets could be very valuable. 
In addition, the use of a single forward-going particle as a probe reduces the uncertainty 
due to jet finding algorithms as well as lowering the minimum angle which can be probed, 
since single particle shower profiles are significantly narrower than comparable energy jets. 
Such considerations led the HI Collaboration to investigate the production of forward vr^s. 
The use of very energetic vrs permits the correspondence between leading particles and 
the struck parton to be used to isolate a region in which BFKL effects could be important. 
In principle any particle species could be used; however, the power of the central tracking 
systems used in the HERA detectors is weakest in this region and the calorimetry in 
general permits both a reasonable identification of vr^s as well as sensitivity to smaller 
angles. 
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Figure 15: Preliminary HI data on production of inclusive jets, found with the 
inclusive kr jet algorithm, in two bins of rj as a function of x. The plots show 
the difference between the data and the NLO QCD prediction. Two predictions 
are compared, one with a normalisation scale = E^, the other with ^\ = Q^. 
The central bands show the effect of hadronisation corrections, while the remaining 
theoretical uncertainty is shown as the outer band. 



The HI collaboration has isolated [|77| a ir^ signal from 5.8 pb~^ of data taken in 1996. 
The 7r° candidates were required to have a transverse momentum in the hadronic CMS 
of greater than 2.5 GeV and a range 2 < < 70 GeV^. They were isolated between 
polar angles of 5° and 25° and required to have energy x^^ > 0.01 ■ Ep, where Ep is the 
proton beam energy of 820 GeV. At such high energies, the two photons from the vr*^ decay 
cannot be separated. Instead, they are identified by a detailed analysis of the longitudinal 
and transverse shape of the energy depositions in order to separate electromagnetic from 
hadronic showers. About 600 7r° candidates were found with > 3.5 GeV. The efficiency 
for detection was around 45%. 



Figure 16 shows the differential cross-section da/dx for priT^^) > 3.5 GeV. Also shown 



are predictions from the RAPGAP and LEPTO Monte Carlo programs, as well as the 
prediction from a modified LO BFKL calculation convoluted with 7r° fragmentation 



functions. The LEPTO model does not give a good description of the data. A 



considerable improvement is given by a model, RAPGAP2.06 |80], which includes resolved 
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Figure 16: The differential cross section as a function of x for forward n^s 
from the HI Collaboration for 2 < < 70 GeV ^. The tt^s were required to 
have pt > 3.5 GeV in the hadronic centre- of-mass system, > 0.01 and angle 
with respect to the proton direction between 5° and 25°. The dashed line shows the 
prediction from RAPGAP and the dotted line that from LEPTO. The solid line 
shows the result of a modified LO BFKL calculation (see text). 

virtual photons in the hard scattering process. Such resolved processes have been shown to 
be important in DIS in some kinematic regions 0] even at moderate Q^. Nevertheless, 
even taking into account the uncertainty caused by varying the renormalisation scale, 
RAPGAP cannot fit the data over the full kinematic range. However, the ARIADNE 
model, which is not shown in the figure, can give a good description of the data, although 
there is considerable arbitrariness in its predictions. The LO BFKL parton calculation is 
in good agreement with the data. However, once again there is a large uncertainty caused 
by a variation in the renormalisation scale of a factor two above and below the nominal 
value, leading to a 60% variation in the prediction. 

Although the agreement of the BFKL model with the data is interesting, overall the 
inherent uncertainties in the various models are such that it is difficult to draw any clear 
conclusion as to the presence of BFKL effects in the data. 



3.3 BFKL tests at the Tevatron 

The Tevatron gives access to a rather different kinematic range in which BFKL effects 
could possibly become important. Here, in the production of high-energy jets, the centre- 
of-mass energy can be much larger than the momentum transfer, Q, so that the jet cross 
section contains large logarithms, ln(s/Q^), which must be summed to all orders. Such a 



summation can be achieved using the BFKL formalism. The DO Collaboration |^ has 



isolated dijet events with very large rapidity separations and measured the cross section 
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as a function of Xi, X2 and Q^, where 1 and 2 label the most-forward and most-backward 
jets, respectively. The longitudinal momentum fractions of the proton and antiproton, xi 
and X2, carried by the two interacting partons are defined as: 

2E 

xi,2 = e^^ cosh(Ar//2) (33) 



where Et^^^Et^) and 771(^72) are the transverse energy and pseudorapidity of the most 
forward (backward) jet, Ar] = r/i — 772 > 0, and f] = {rji + ?72)/2. Thus, by measuring the 
cross section at the largest accessible values of Ar/, the separation in x of the colliding 
partons is maximised, thereby optimising the phase space for gluon ladders with strong 
ordering in x, i.e. BFKL evolution, to be produced. The momentum transfer during the 
hard scattering is defined as: 

Q = a/ Et^Et2 
The BFKL prediction [^] for the cross section is: 



"Tbfkl oc — — (34) 

where aBFKL is the BFKL intercept that governs the strength of the growth of the gluon 
distribution at small x, which in leading logarithmic approximation (LLA), is given by: 

. «5(Q) 12 In 2 

01 BFKL - 1 = (35) 

TT 



The cross section is a convolution of the probability to find the interacting partons with 
given values of x inside the proton, together with the partonic scattering cross sections 
that contain the BFKL effects of interest. In order to remove the dependence of the cross 
section on the parton distribution functions, it is convenient to measure the cross sections 
at different centre-of-mass energy and take the ratio at fixed xi,X2 and Q^] this also 
reduces some correlated systematic effects stemming from the detector characteristics. 



The ratio of cross sections can be obtained from equation [34 



= / , X = 7 = , = (3o) 



Clearly, the greater the difference between the two centre-of-mass energies, the greater 
the effective variation in Arj between the two data sets and thereby the larger the BFKL 
effects should become. 

The data samples used in the DO analysis was taken at a/s = 1800 and 630 GeV. The 
Xi,X2 ranges were restricted to lie between 0.06 and 0.22 in order to avoid detector biases. 
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Figure 17: T/ie raizo o/jei cross sections with rapidity gaps between them mea- 
sured at ^Js = 1800 and 630 GeV by the DO Collaboration. The data points are 
plotted as filled circles for rapidity gaps greater than 1, and 2, as a function of the 
mean rapidity gap at 630 GeV. The light shaded circles represent the predictions 
from HERWIG while the square represents a calculation based on the LLA BFKL 
approach. The thick error bars denote the statistical uncertainty, while the thin 
error bar denotes the statistic and systematic uncertainties added in quadrature. 



Only one bin was used, with 400 < < 1000 GeV^. Figure [T^ shows the ratio of the 
cross sections at the two centre-of-mass energies as a function of the mean separation in 
r] in the 630 GeV data set. Also shown on the plot are the predictions from LO GLAP 
evolution, from the HERWIG Monte Carlo, and from the BFKL calculation. The data 
agree with none of the models, showing a much steeper increase than predicted even by 
the BFKL calculation. However, the steep rise is certainly more in accord with the BFKL 
prediction than that of LO GLAP evolution; what the rise in the HERWIG predictions 
means is unclear, at least to this author. In conclusion, the situation is yet again confused. 

To summarise this section, despite having examined specific exclusive processes in which 
the effects of BFKL evolution are expected to be maximal, the current experimental sit- 
uation is that there is little firmer evidence for deviations from standard GLAP evolution 
than there was in the inclusive F2 data. 



4 Interpretation and Models 

While the standard perturbative QCD GLAP evolution reigns supreme at medium x and 
high Q^, the area of low x is a particularly rich and complex area in which this perturbative 
QCD anzatz meets and competes with a large variety of other approaches, some based 
on QCD, others either on older paradigms such as Regge theory or essentially ad-hoc 
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phenomenological models. In this section a far-from-exhaustive survey of such models 
is undertaken. Often it is convenient to concentrate on a particular model in order to 
confront specific predictions with the data. This in no way necessarily implies that such 
an example model is the best one available, or that others can be ignored in its favour - 
rather it is an attempt to illustrate generic characteristics of particular classes of approach 
to the elucidation of the data. 

In this section the older, Regge-based approach to the understanding of F2 is examined, 
followed by models that give simple parameterizations of the structure functions based, 
more or less loosely, on pQCD. A brief overview of the current state of global fits to the 
parton distribution functions is then carried out, followed by a discussion of new data from 
ZEUS on the logarithmic derivatives of F2 and possible implications for QCD evolution. 
Dipole models, in particular that of Golec-Biernat & Wiisthoff p4|,p5| are discussed in 
some detail and compared to the data. 



4.1 Regge-based models of F2 



Regge theory [[72| , p6| , [87[1 has a long and distinguished history in particle physics as a con- 
ceptual basis linking bound-state spectra, forces between particles and the cross-section 
behaviour as a function of energy over a wide variety of processes through the analytical 
properties of high-energy scattering amplitudes. The basis of the theory is that solutions 
of the Schrodinger equation for non-relativistic potential scattering can be solved in terms 
of a complex angular momentum variable, j. For many potentials, the only singularities 
of the scattering amplitude are poles in the complex angular momentum plane, known as 
'Regge poles'. Somewhat surprisingly, such non-relativistic concepts have proven to be 
very useful in particle physics. To take the specific example relevant to this discussion, 
any total cross section having a power-law dependence on the centre-of-mass energy lends 
itself to a simple explanation in terms of Regge poles corresponding to the exchange of 
specific particles. The most important poles are those corresponding to the exchange of 
vector and tensor mesons, which give a pole at j = 1/2, and that near to j = 1, which 
corresponds to a particle with the quantum numbers of the vacuum that has never been 
observed as a final-state particle. This particle is known as the Pomeron, and is thought 
of as the exchanged particle that mediates elastic and diffractive scattering. In fact, if the 
pole is assumed to be at j = 1.08, known as the 'soft Pomeron' singularity, then it gives 
a good description of the cross sections of a wide variety of hadronic processes |^ . 

The advent of HERA, allowing access to large in both soft and hard processes, offered a 
fertile ground for the application of Regge theory, since it works par excellence at very high 
energies; for example, in DIS when W'^ is much greater than any other invariant. Indeed, 
equations O and E3| make it clear that, at low x, F2 can be expressed in terms of high- 
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energy transverse and longitudinal cross sections. It was therefore interesting to discover 
that the strong rise in F2 at low x discovered at HERA, which can be parameterised as 
W^, was incompatible with the soft Pomeron intercept that had worked so well in 



many other reactions. Donnachie and Landshoff suggested that a good description 
could be achieved if a further simple Regge pole were assumed at j = 1.435, the so-called 
'hard Pomeron'. They fit to the form: 



]X 



(37) 



and, as can be seen from figure [T§, such an ansatz does indeed give a good description of 
the F2 data at low x (and indeed at higher x as well). 



16-05 




Figure 18: The HI published data on F2 plotted divided by in different bins 
as a function of x. The lines show Donnachie and Landshoff 's fits using Reggeon, 
soft and hard Pomeron contributions. 



Donnachie and Landshoff also looked at the ZEUS data on the charm structure function, 
^2'^ 0]' Here, presumably because the mass of the charm quark provides a hard scale, 
there is no requirement for the vector and tensor meson or the 'soft Pomeron' pole, leaving 



only one term in equation 3/. The fit shown in figure 19 to the 'hard Pomeron' term only 
gives a good fit to the data. 

Despite the successes of this approach, there are several major drawbacks. Firstly, the 
fi{Q^) in equation 3? cannot be predicted from Regge theory. Having added one ex- 
tra Pomeron, there is no logical reason not to add others when HERA data from other 
channels plj-p^ and improved accuracy imply deviations from a 'universal' two-Pomeron 
ansatz. Indeed, more complex singularities than poles can also be added Il9^-p6i leading to 
behaviour more complex than simple powers of W. Thus, the Regge description can have 
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Figure 19: ZEUS data on plotted divided by as a function of x. The lines 
show the Donnachie-Landshoff prediction with only the hard Pomeron term in the 
fit. 

many essentially arbitrary parameters and risks becoming merely a phenomenological pa- 
rameterisation. The great range of seemingly disparate phenomena that can be explained 
in the Regge framework nevertheless implies that any more comprehensive theory of the 
strong interaction, viz. QCD, must somehow 'explain' or assimilate Regge concepts in a 
natural way. This line of investigation is a fruitful one currently being actively pursued 
and also leads naturally into the next section. 



4.2 QCD-inspired parameterisations 



The double-logarithmic limit of QCD, in which both and 1/x become very large, has 
long been known to imply |^ that F2 should fit to the form: 



F2{x, oc exp{v/(48//3o)ln(l/x)lnlng2} 



(3^ 



where /3o is the standard Renormalisation Group Equation (5 function. Ball and Forte 
p9|] showed that indeed the HERA data were well represented by this form, which they 
refer to as 'double-asymptotic scaling'. 

Other authors have examined functional forms that are related to the double-asymptotic 
scaling of equation Buchmiiller and Haidt | |10CI| | set out the theoretical region of 
validity of the so-called 'double-logarithmic' scaling regime. The functional formQis: 



F2(x, Q"^) = m In 



Q 



ln(^ 

X 



(39) 



^ Modified [101| from the original form [102| proportional to InQ^/Qp in order to take account of new 
data at very low Q^. 
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where m is a constant to be determined from the data. Figure ^ shows [ 103 | HERA dataQ 
plotted as a function of to this functional form. A linear fit | |103|| is remarkably good. 




Figure 20: ZEUS and HI F2 data plotted against the scaling variable ^ (see 



equation p^ j 



This functional form is appropriate to the Regge hmit of fixed and x — 0, and indeed 
it corresponds to the first term of the more general solution to the double-asymptotic 



form viz. a summation of all terms of the form | 100 | 



1 Q% 1 



(40) 



which rises faster than any power of Inl/x as x falls. In the more general case in which 
both and 1/x become large, or at sufficiently small x, the double-asymptotic form 
ought to become more appropriate. In fact, somewhat surprisingly, it would seem that 
even at very small x the double-logarithmic form gives a good fit to the data. 

A parameterisation clearly related to both the double-asymptotic and double-logarithmic 
ones has recently been used by Erdmann ||104|| , in particular to facilitate comparison 
between the proton, photon and Pomeron structure functions. It has the form: 



F2(x,g2) = a{x) 



In 



A2 



(41) 



where A is the QCD scale parameter. Figure ^ shows that this form gives a good fit to 
the preliminary HI F2 data. 

It is also true for the published ZEUS and the BCDMS data, even at high x. In some 



sense equation 41 represents the alternative approximation to the Regge limit of Haidt, 



The data at the lowest was stiU prehminary at the time of this conference. 
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Figure 21: Preliminary HI F2 data fit to the form of equation ^1 



i.e. fixed x and large. The fact tliat botli a and k are functions of x liowever allows the 
In 1/x-like terms in the sum of equation 40 to be approximately taken into account, so that 
the parameterisation also works well at low and x. It can be seen by comparison with 
equation |T8| that a{x) is proportional to the charged parton distribution at a particular 
value of Q^, viz. that for which InQ^/A^ = 1, and that k is related to the scale-breaking 
of F2 and therefore to gluon radiation from the partons. 

Figure ^ shows the values of a and k as a function of x. The valence quark distribution 
at high X can be very clearly seen, as well as that fact that for the low-x regime of interest 
here, a{x) is approximately constant. This implies that the scale-breaking represented 
by K is what drives the increase in F^. Figure shows that k increases more or less 
linearly from the negative scale-breaking at high x until at least x^lO"^, at which point 
it seems to level off. The implications of this observation will be discussed in more detail 
in section |4.4| in conjunction with the ZEUS data at very low x. 



4.3 Global fits 



The extension of the kinematic range and the high-precision data on F2 from HERA 
provided a substantial impetus to the determination of parton distribution functions 
via global fits to a wide variety of data. The major approaches are due to the CTEQ 
group |05l, Gliick, Reya and Vogt (GRV) |lOi and Martin et al. (MRST) |10^. In gen- 
eral all groups fit to data from fixed target muon and neutrino deep inelastic scattering 
data, the HERA DIS data from HERMES, HI and ZEUS, the W^-asymmetry data from 
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Figure 22: The values of the coefficients in equation^^ as a function of x after 
fits of this functional form to BCDMS, HI and ZEUS F2 data (filled symbols). 
The open symbols correspond to photoproduction and diffractive data that are not 
discussed in the text. 

the Tevatron as well as to selected process varying from group to group such as prompt 
photon data from Fermilab as well as high-i?r jet production at the Tevatron. The dif- 
ferent data sets give different sensitivity to the proton distributions depending on the 
kinematic range, but together constrain them across almost the whole kinematic plane, 
with the possible exception of the very largest values of x, where significant uncertainties 
still remain ||108 |. 



The approach of GRV is somewhat different from that of the other two groups. They 
utilise the fact that, as ^0, parton distributions are fully constrained by the charge 
and momentum sum rules. By assuming valence-like distributions for the quarks at a 
very low starting Q^, in principle the gluon and sea distributions can be generated purely 
dynamically. However, it is found that such a procedure generates parton distributions 
which are too steep as x decreases. Instead they input 'valence-like' distributions for both 
quarks and gluons fixed by high-a; data at a larger though still very small Q^. The starting 
value, Qq, is determined by the point at which the input gluon distribution is of the same 
order as the input u valence quark distribution and is ~ 0.5 GeV^ in NLO QCD [ 109 



Although there are quite large uncertainties on the value of Ql and on the valence-like 



distributions assumed at Ql, the effect of these is suppressed in the comparison with the 
high-Q^ HERA data by the long evolution distance. In general, the GRV parameterisation 
gives good fits to the HERA data, as shown in figure although as — > Qq the fit 
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becomes worse, as would be expected from the formalism. In addition, however, GRV at 
NLO has difficulties in fitting the logarithmic derivatives of F2 for values of a; < 10^^ ||11CI|| 
(although see section [4.4] ). 
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Figure 23: The GRV98 fit to the ZEUS and HI F2 data in the low-Q"^ region in 
bins of plotted as a function of x. 

The approaches of CTEQ and MRST are basically similar, although they differ both 
in the data sets used as well as in the fitting procedure and the technical details of the 
theoretical tools used, e.g. the treatment of heavy quarks in DIS. In their latest fits CTEQ 
prefer to omit the prompt photon data because of the uncertainties in scale dependence 
and the appropriate value for the intrinsic required to fit the data. Instead they use 
single-jet inclusive distributions to constrain the gluon distribution at large x. In 
contrast, until their most recent publication, MRST retained the prompt photon data, 
giving alternative PDFs depending on the value for the prompt-photon intrinsic used. 
Both groups parameterise the parton distributions in terms of powers of x and (1 — x) 
leading to fits with many free parameters. The MRST NLO parameterisation of the gluon 
is shown below as an example: 

xg = AgX-^' (1 - (1 + eg^/^ + -fgx) (42) 

where Ag, Xg,rig, eg and 7^ are free parameters in the fit. The treatment of the d/u 
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ratio at high x has recently been addressed by Yang and Bodek | 108 |, who point out 
that deuterium binding corrections should be applied to the NMC F^jF"^ data. Such 
corrections give good fits in the global analyses, except to the uncorrected NMC data 
themselves. The PDFs determined from the CTEQ5M fit are shown in figure ^ at = 
25 GeV^. The steep rise of the gluon and sea distributions as x falls is evident. 




Figure 24: The PDFs resulting from the CTEQ5M fits at = 25 GeV^. 

A long-standing problem with the various global PDFs has been the fact that no error 
was associated with the central values. The difficulties associated with producing such 
errors from a multi-parameter fit to many data sets with differing correlations are cer- 
tainly formidable. It is therefore an extremely important and welcome development that 
Botje has recently produced for the first time PDFs with associated error matrices [ p. 11 



Figure ^ shows the valence- and sea-quark distributions together with that of the gluon 
from Botje's analysis. 

The fit utilises a more restricted range of data than the CTEQ and MRST fits, using 
the HI and ZEUS F2 data together with the fixed target muon and neutrino data; the 



Drell-Yan data from E866 ||112|| are used to constrain the u — d distribution. Despite the 
more restricted data sets used, the results of the fit are very compatible with the most 
recent fits of CTEQ and MRST. The importance of the error matrices produced by this 



fit can be illustrated by the example shown in figure Here the effect of the errors on 
the uncertainty in the prediction of various cross sections and ratios in the ZEUS DIS 
dijet analysis ||113|| is shown. The very large difference in the estimated error, depending 



on whether the correlations between the parameters in the PDFs are taken into account 
of not, is striking and of the highest importance in a realistic calculation of the error on 
quantities such as a^. 

All of the parameterisations discussed above were carried out in the framework of NLO 
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Figure 25: The valence quark, sea quark and gluon PDFs resulting from the Botje 
fit. The bands show the uncertainty associated with each PDF. 
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Figure 26: The uncertainty of three quantities used in the ZEUS determination 
of Os from dijets in DIS. The cross section differential in , the dijet differential 
cross section, and the ratio of the dijet to the total cross section are shown. The 
shaded area shows the effect of taking the correlated errors produced by the Botje 
fit properly into account, the hatched area that of ignoring the correlations. 

QCD. With the increasing precision of the DIS data, as well as the need for accurate pre- 
dictions of cross sections at the LHC, the need for next-to-next-to-leading-order (NNLO) 
fits is obvious. The first steps in this regard have already begun, and some moments of the 
NNLO splitting functions have already been calculated [p.l4|| . Using this with other avail- 



able information, van Neerven and Vogt ||115| , |116|] have produced analytical expressions 
for the splitting functions which represent the slowest and faster evolution consistent with 
the currently available information. The MRST group has recently used this information 
to investigate NNLO fits to the available data | 117 |. Such an analysis requires some 
changes to the parameterisations used, so that for example the NLO parameterisation of 
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the gluon of equation 02] becomes: 



xg{x,Ql) = AgX (l+Egy/x + ^gx) 



A'g x-^'^ (1 



xf'^ 



(43) 



primarily in order to facilitate a negative gluon density at low x and low Q^, which, 
although conceptually somewhat bizarre, is nevertheless preferred by the fits, even at 
NLO. The results of the 'central' fit, between the extremes of the van Neerven-Vogt 
parameterisation, is shown in figure p7. 




Figure 27: The MRST 'central' NNLO fit to DIS data. The solid Ime shows 
the NNLO fit, while the NLO fit is shown by the dashed line and the LO fit by the 
dotted line. The data are from HI, ZEUS and the fixed target experiments and are 
plotted in x bins as a function of with an additive constant added to the data of 
each X bin to improve visibility. 



There are also changes of the LO and NLO fits with respect to earlier publications, in 
as much as MRST now follow CTEQ in using the Tevatron high-E'T- data rather than 
the prompt-photon data, and preliminary HERA F2 data has been included in the fit. 
There is a marked improvement in the quality of the fit in the progression LO NLO 
— > NNLO, in particular in terms of the NMC data. The size of higher-twist contributions 
at low X also decreases, so that at NNLO is it essentially negligible. The effect of going 
to NNLO on the PDFs themselves is highly non-trivial. This is illustrated in figure E3, 
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where the quite major changes in F^, particularly at low x, are evident. There is also 
a large variation depending on the choices made in the parameter space allowed by the 
partial NNLO ansatz. Indeed, the GLAP approach is not convergent for < 5GeV^, 
which may well be due to the neglect of important In 1/x contributions. However, the 
instability seen at low soon vanishes at higher Q^. 




Figure 28: The Fl structure function from the MRST fits, taking into account part 
of the NNLO corrections in four bins of Q"^ as a function ofx. The solid line shows 
the 'average' of the parameter space available to choose the NNLO parameters, while 
the dashed- dotted lines show the two extreme possibilities. The NLO fit is indicated 
by the dashed line while the LO fit is indicated by the dotted line. 



Thorne has indeed investigated the question of incorporating In 1/a: terms in the splitting 
functions by incorporating the solution of the NLO BFKL kernel using a running coupling 
constant [|118|| . The results are shown in figure |29|. It is clear that the inclusion of 



the BFKL terms does indeed give an improved fit compared to the 'central' NNLO fit, 
particularly at the lowest and x. This may be one of the first unambiguous indications 
of the importance of BFKL evolution; if so, it is rather surprising that it has occurred 
in the analysis of the inclusive data, rather than the exclusive channels, expected to be 
more sensitive, that were examined in section |[ 

In conclusion, there have been major advances in the field of parton distribution func- 
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MRST NNLO and NLO fits , x = 0.00005 - 0.00032 
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Figure 29: The NNLO MRST fit, modified by Thorne to include a NLO BFKL 
kernel. The 'standard' NNLO fit is shown as the full line, the dashed line shows 
the NLO fit and the dotted line the Thorne BFKL modification. 

tions and global fits in the last year. Not only is there now a parameterisation which 
gives produces associated error matrices, which is of first importance in the treatment 
and extraction of experimental results, but also the first attempts to incorporate NNLO 
corrections into the fitting has begun. In the latter case, it is clear that there is still a great 
deal of work required before there is a real understanding of the effects of a full NNLO 
treatment; not the least of the work is in the onerous task of deriving all the necessary 
NNLO terms. It may still be premature [|120|| to worry too much about the somewhat 



strange behaviour of the NNLO gluon density and Fl, until a full NNLO treatment is 
possible. Nevertheless, the increased precision of the data becoming available and the 
rapid theoretical developments combine to make the subject of global PDF fitting and 
structure functions both topical and interesting. 



4.4 F2 and its derivatives 

With the publication of the final data from the very \ow-{Q^, x) region measured with the 
Beam Pipe Tracker (BPT) as well as the latest high-precision F2 data, ZEUS now has 
precise data over a remarkable six orders of magnitude in x and Q^. This data is shown 
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in X bins as a function of InQ^ in figure 30, togetlier with fixed target data from NMC 
and E665, wliicli extends tlie range in tlie direction of medium x and Q^. 




Figure 30: Compilation of ZEUS F2 data, both published and preliminary from 
the 1996-97 data sample in x bins as a function of . Each x bin is shifted by an 
additive constant for ease of visibility. Data from NMC and E665 is also shown. 
The dotted lines show lines of constant W , while the solid lines are fits to the form 



of equation ^ . 



The availabihty of this very wide range of precise data makes possible quahtatively new 
investigations of models that describe F2. As discussed in section pl2| , the logarithmic 
derivative of F2 is directly proportional to the gluon density, which in turn is by far the 
dominant parton density at the small values oix of interest here. It is therefore interesting 
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to examine the behaviour of such logarithmic derivatives as a function of both x and Q^. 
Plots of dF2/d\nQ'^ as a function of x were first presented for the ZEUS data by Caldwell 
at the DESY Theory Workshop in 1997 and subsequently published by ZEUS [ p.21|| , and 



led to much comment in the literature. The range and quality of the data available at 
the time meant that severe restrictions were placed on how the data could be binned and 
parameterised. These restrictions led to several erroneous suggestions that the features of 
this plot were a consequence of trivial kinematics. The quality and range of the currently 
available data now permits a much better-defined procedure to be followed in constructing 
plots of the logarithmic derivative. 

The data shown in figure [3^, particularly in the lower-x bins, are clearly inconsistent 
with a linear dependence on InQ^, as was pointed out for the preliminary HI data by 



Klein |32|. The solid curves on the figure correspond to fits to a polynomial in InQ of 



the form: 



F, = A{x) + B{x) (logio Q') + C{x) (logio Q')' (44) 



which gives a good fit to the data through the entire kinematic range. The dotted lines 
on figure ^ are lines of constant W. The curious 'bulging' shape of these contours of 
constant W in the small-x region immediately implies that something interesting is going 
on there. Indeed, simple inspection of figure |3^ shows that the slope of F2 at constant 
W begins flat in the scaling region, increases markedly as the gluon grows and drives the 
evolution of F2 and then flattens off again at the lowest x. 



This behaviour is made clear and explicit in figure ^ which shows the logarithmic deriva- 
tive evaluated at (x, Q^) points along the contours of fixed W shown on figure |30| according 



to the derivative of equation 44, viz.: 



g^ = BW + 2CWl„g,Q^ (45) 

where the data are plotted separately as functions of InQ^ and Inx. The data was 
adjusted where necessary to the appropriate Q^, x, W-hin by using the ALLM parameter- 



isation |122|. The error bars on the points are evaluated from the errors on the parameters 



in equation ^ and consist of the statistical and systematic uncertainties added in quadra- 
ture. The correlations on the errors are, however, not taken into account, so that the error 
bars shown are slight over-estimates. 

The turn-over in the derivatives in all W bins is marked, and confirms the similar feature 
seen in the original ZEUS plot, but now with much better defined kinematic conditions. 
When plotted as a function of InQ^, the maximum in the derivative moves to larger as 
W increases, while as a function of In x, the maximum moves to smaller x asW increases. 
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Figure 31: The logarithmic derivative of the ZEUS F2 data in six bins of W, 
plotted as a function of Q"^ and x. 

It is of interest to speculate at this point as to what dynamical mechanism might be caus- 
ing the behaviour exhibited in figure |31|. Since the logarithmic derivative is proportional 
at leading order to the gluon density, the obvious inference that can be drawn from the 
data is that the rise in the gluon density at low x begins to soften and eventually to fall 
as X decreases. Indeed, several indications of such an effect have been discussed in earlier 
sections of this talk. Such an effect is by no means, as will be seen below, necessarily an 
indication of deviations from GLAP evolution. Nevertheless, such a fall in the gluon den- 
sity as X falls is a natural consequence of many models of parton saturation or shadowing, 
so that it is of interest to explore their features in more detail at this point. Before begin- 
ning however, it is important to emphasise that the relative emphasis on dipole models 
in this talk is not an indication that they are necessarily 'correct', or even necessarily 
give a better description of the data than other models, such as the standard twist-two 
QCD descriptions. Neverthless, they do have several attractive features, in particular the 
rather natural way in which they can lead to a unified description of diffraction and deep 
inelastic scattering, which makes it useful to discuss their features in some detail here; 
not least since in general their concepts are less familiar to the average particle physicist. 
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4.4.1 Dipole models and shadowing 



Dipole models of DIS have a long history. The basic idea is to transform the 'normal' 
way of looking at DIS, which considers a virtual photon to be emitted from the incoming 
lepton and to collide with a parton emanating from the proton, by transforming to a 
topologically equivalent process in which the virtual photon splits into a quark-antiquark 
pair. These two descriptions are related by a Lorentz transform, since the 'normal' view 
of partons evolving inside the proton is appropriate to a frame such as the Breit frame 
or the infinite-momentum frame, whereas the dipole picture is more appropriate to the 
rest frame of proton. In the rest frame of the proton, the virtual photon splits into a 
quark-antiquark pair, or dipole, well downstream of the proton. The formation time of 
the dipole in the proton rest frame is related to the uncertainty in the energy of the pair 
by Tqq ~ 1/AE, which, in the limit of small x becomes [|123|] Tgg ~ l/{xMp), where Mp 



is the proton rest mass. Since the distance between the formation of the dipole and the 
interaction with the proton implied by this lifetime is much larger than the proton radius, 
the transverse size of the dipole can be considered fixed during the interaction. Thus, 
for small x, the deep inelastic process can be considered semi-classically as the coherent 
interaction of the dipole with the stationary colour field of the proton a long time after 
the formation of the dipole. In such a frame the dipole does not evolve a complex parton 
structure, which is considered to take place inside the proton. 

It is clear that the formulation of DIS in this dipole picture provides a direct link between 
the processes of deep inelastic scattering and diffraction. The fully inclusive structure 
functions sum over all possible exchanges between the dipole and the proton, dominantly 
one- and two-gluon exchange, whereas diffraction is produced by the exchange of 2 glu- 
ons in a colour-singlet state. This deep connection between these two processes leads 
to non-trivial predictions which do indeed seem to be borne out by the data. They 
have been investigated by several authors, including Golec-Biernat and Wiisthoff ||8^ and 
Buchmiiller, Gehrmann and Hebecker ||50| . 

Qualitatively, the interaction of the dipole with the colour field of the proton will clearly 
depend on the size of the dipole, which is proportional to Q~^. If the separation of the 
quark and antiquark is very small, the colour field of the dipole will be effectively screened 
and the proton will be essentially 'transparent' to the dipole. At large dipole sizes, the 
colour field of the dipole is large and it interacts strongly with the target and is sensitive 
both to its structure and size. Such considerations lead naturally to some qualitative 
understanding of the process of deep inelastic scattering and saturation illustrated in a 
simple one- dimensional model in figure |32|. Here a one- dimensional distribution of partons 
inside the proton is considered in two limiting cases. In the first, labelled as 'scaling', the 
typical size of the probing dipole is much smaller than the mean separation of the partons. 
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Figure 32: Schematic view of the 'scaling' and 'saturation' regimes in DIS. The 
relative sizes of the dipole (proportional to 1/Q) and the mean separation between 
partons control the behaviour. 

Rq, so that the probabihty of interaction is given by the ratio between the mean size of 
the dipole and the mean separation of the partons, i.e. 1/{QRq). The cross section is 
thus proportional to l/((5-Ro)^; so that the structure function is independent of Q^. In 
the other case, labelled as 'saturation', the size of the dipole is large compared to the 
mean separation of the partons, in which case the size which determines the interaction 
probability is simply the size of the probe. Thus, for a given Q, the cross section 'saturates' 
to a constant value. More generally, when the parton density is such that the proton 
becomes 'black' and the interaction probability is unity, the dipole cross section saturates 
for all Q and hence the structure function becomes proportional to Q^. In the Breit-frame- 
like picture, this is equivalent to a situation in which the individual partons become so close 
that they have a significant probability of interacting with each other before interaction 
with the probe. In the case of gluons, such interactions lead to two one branchings and 
hence a reduction in the gluon density. Such a picture was the basis for many of the early 
developments in this area, in particular the formulation of the modified GLAP evolution 
equation including absorptive effects by Gribov, Levin and Riskin and Mueller and 
Qiu [EO] , as embodied in the GLR equation. 



In order to look somewhat more quantitatively at the implications of these ideas, it is 
necessary to specialise to a particular model. The dipole model of Golec-Biernat and 
Wiisthoff [§5|S^ (G-B&W) is selected to be discussed in detail. This does not of course 
imply that many other models P2|- pn| , p!^^l 2 71 | both older and more recent, are not equally 
capable of describing the data. In particular, most dipole models share many of the 
characteristics of the G-B&W model, at least at the rather broad-brush level appropriate 
to this discussion. 

The interaction between a dipole with a definite transverse separation ^ at a fixed impact 



parameter b and the proton can be considered very generally |52] in terms of an S'-matrix 
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element ^(.^jb). The form assumed for this ^-matrix element contains the basic physics 
of the model in question. In the G-B&W model, it is assumed that the impact-parameter 
dependence can be factorised and integrated over and that the remaining dependence 
on the separation can be approximated by a Gaussian. Explicitly, the cross section for 
transverse and longitudinal photon is given by: 

1 

aT,L(x,Q') = Jd^^J dz\<i/T,L{z,r)\^a{x,r^) (46) 



where '^t,l are the light-cone wave-functions for the photon, which are functions of the 
fractional momentum of the virtual photon taken by quark, z, and the separation of the 
quark and antiquark r. The photon wave-function has the form: 

1^^ = S E 4 + (1 - ^'Klier) + m)Kl{er)] (47) 



\'^l{z,-)\" = [4Q^-^(1 - ^fKl{er)] (48) 



where Kq and Ki are McDonald functions and 

= z{l-z)Q'^ + m) (49) 

where is the mass of the quark in the dipole. Neglecting for the moment the fermion 
mass, the fact that the argument of Kq and Ki is er implies that the 'effective' size 
of the dipole configuration is proportional to 1 / {Q\/ z{l — z)}. Thus, the fact that the 
longitudinal wave-function from equation ^ is proportional to z{l — z), whereas the 
transverse wave- function in equation ^ is proportional to 2;^ + (1 — 2;)^ implies that 
the larger configurations, when 2 or (1 — 2) — > 0, are suppressed for the longitudinal 
photons ||128|| . For large dipole configurations, the integral over z in equation ^ picks 



up contributions only from the end-points, in which either the quark or the antiquark 
carries essentially all of the photon momentum; such configurations are therefore known 
as 'aligned'. Since the colour field and hence the interaction probability is lower for 
smaller dipoles, the dipole cross section is dominated in most areas of phase space by the 
transversely polarised component of the virtual photon. 

The sub-process cross section, a, in equation ^ is related to the 5'-matrix element dis- 
cussed above, and is assumed in the G-B&W to have the form: 

a(a;,0 = a,g{f^) (50) 
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where: 

9 (r^) = 1 - e~'' 



(51) 



2i?o(x) 

and: 



(52) 



Ro{x) = -^{-] (53) 

These definitions contain the essential dynamics of the G-B&W model. At large 'rescaled' 
dipole sizes, f, (7 — > constant and the cross section saturates. For small f, the cross-section 
increases quadratically with f, which, from equations ^ and |5^, implies an rise as 
seen in the data. In order to pick up the dependence it is necessary to do the integral 
in equation ^ for the transverse component. At small er, the McDonald functions can 
be approximated by: 

Ko{er) ~ In— (54) 
er 

i^i(er) ~ — (55) 
er 

while at large er they are exponentially suppressed. Thus, it is clear that the dominant 
contribution to the integral comes from the Ki term for er < 1. This corresponds to the 
small f case discussed above so that for 'small' dipoles, i.e. r < 1/Q « Rq, for which 
er < 1 is automatically satisfied, the saturation radius becomes: 



Substituting in equation E6| using equation the integral collapses to: 







since the z integral can be factored out since the e terms cancel and the r integral is limited 
to an upper limit of 1/Q by construction. At constant x, therefore, F2 (oc Q'^a^'^) exhibits 
scaling. By analogy it is easy to see that for 'small' dipoles in which the characteristic 



size 1/Q > Rq, the integral in equation |56| must be split into two parts, in which a is 



quadratic in r for small values of r and constant for large r, i.e. 

ar ~ / ( — j ctq^ + / c/r^ ( — ) (Tq oc ao + ao In ( „^ ) (57) 



Rl J V^V \Q^Rl 



Rl 
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which predicts that F2 is proportional to (modified by a slow logarithmic dependence). 
Thus it can be seen that the presence of an additional length scale in the problem, the 
saturation radius, leads to the prediction that for sufficiently large dipoles (i.e. small Q^), 
F2 will become proportional to Q^. The boundary between these two types of behaviour 
is the so-called 'critical line' (given by 1/Q = Rq{x)), which clearly depends on x. With 
increasing W, the transition occurs for smaller x and larger Q^. 

Having given a broad-brush overview of the main implications of the G-B&W model, 
it's predictions for the logarithmic slope of F2 can be investigated. The more detailed 
treatment in Golec-Biernat and Wiisthoff modifies the conclusions of equations ^ 



and ^ by the inclusion, among other factors, of 'large' dipole pairs together with the 
longitudinal contribution, to give: 



a''''{x,Q' 




x_y 

x'qJ Ql 



In 



g2 ^ 



(58) 



where the primes denote that the constants are to be optimised by a fit to the available 
data. The salient characteristics of equations ^ and ^ remain, although equation ^ has 
acquired a logarithmic modification. The first term therefore governs the behaviour at 
high Q^, while the second term is dominant at low Q^. Multiplying equation ^ by to 
convert it to F2 and taking the logarithmic derivative leads to: 



dF2 _v 
~ X 



for high Q^(>> l/i?o) ^'^ derivative acquiring a term proportional to 

for low (5^(~ V-^o) fixed x, thereby reducing the size of the derivative. Thus the 
expected power-law growth at low x is seen for high Q^, where the logarithmic derivative 
in LO GLAP treatment is proportional to the gluon density, while at small the leading 
behaviour of both the derivative and F2 becomes proportional to Q^. This implies a 
maximum in the logarithmic derivative, as seen in the data of figure Qualitatively, 
therefore, the G-B&W model can describe the ZEUS data, as shown by figure ^ which 
contains the curves from the original Golec-Biernat and Wiisthoff publication, which 
indeed is qualitatively in agreement with the ZEUS data. In particular, the movement of 
the maximum with and x as W changes is quite well reproduced, but there are clear 
differences, particularly at higher and lower W. These are in regions in which the 
model has known problems, and it could well be that a fit to the ZEUS data, which were 
not available at the time of the original paper, would improve the agreement. 

It is also of interest to examine the logarithmic derivative at fixed rather than fixed 
W as a. function of x. This is shown for the ZEUS data in figure 53. The derivative 
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Figure 33: Curves showing the G-B&W model predictions for the logarithmic 
derivatives of F2 in bins of constant W as a function of and x, compared to the 
ZEUS data. The curves are only plotted for x < 10~^, the limit of validity of the 
model. 

is relatively straight as a function of In x, and exhibits a slow change with at larger 
Q^, which becomes rapid for < 2 GeV^. This behaviour is simply a reflection of the 
'valence-like' gluon behaviour at low in which the gluon density in NLO QCD fits 
falls rapidly to zero while the sea remains non-zero. It is also reproduced by the G-B&W 
model. The logarithmic derivative of equation multiplied by at fixed shows 
that the slope changes from being proportional to x""^^' to Inx"'*'', so that there is simply 
a change in slope rather than a turn-over. Moreover, the G-B&W 'critical line', which 
predicts the position of the transition to saturation behaviour, is much steeper in x than 
in Q^, so that the transition point for fixed is generally at an x outside the kinematic 
region of the data. The exception is the = 0.75 GeV^ data, where the transition is 
predicted to occur at around x ~ 5 ■ 10"'^. Such a change in slope can certainly not be 
ruled out by the data. 

Although the qualitative agreement of the G-B&W saturation model with the ZEUS data 
is intriguing, there are many other possible explanations. Several other saturation and/or 
dipole models can describe the general trends of the data [ 128 , 129 . It was already re- 



marked that the parameterisation of Haidt of equation 39 also gives rise to a turnover in 
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Figure 34: The logarithmic derivatives of the F2 data in bins of constant as 
a function of x. 

the logarithmic derivative. Furthermore, several NLO QCD analyses seem able to repro- 
duce the turn-over and other general features of the data. Bliimlein ||130|| has used the 
GRV framework to fit qualitatively the ZEUS data. Roberts ||131|| has produced logarith- 



mic derivatives using the MRST fits; although they do produce a turn-over, its position 
and the lower and x slopes do not agree well with the data. This is scarcely surprising 
since the whole parton picture must already be questionable in such a kinematic region. 



and there may well be important higher-twist contributions. Nevertheless, Thorne [02 



including the BFKL-motivated modification of the splitting functions discussed in *sec- 



tion 44.3, section 4.3, has produced modified MRST fits that give an improved fit to the 



ZEUS data. Both results are shown in figure 35 



From the above it is clear that the current ZEUS data can certainly not be used to claim 
evidence for saturation effects in the HERA kinematic range. Higher-precision data will 
certainly help in distinguishing competing explanations. However, it does not seem likely 
that the existence or otherwise of parton saturation can be unambiguously established 
at HERA, at least from studies of the logarithmic derivatives of F2 alone. One problem 
is that the centre-of-mass energy of HERA means that the interesting areas at low x 
in which the saturation effects become large is necessarily at < 5 — 10 GeV^, where 
complications from higher-twist effects are inevitable and indeed the whole parton picture 
at some point ought to break down. The only way to improve the situation would be to 
move to a higher energy machine - for instance the proposed THERA option of colliding 



TESLA and HERA [|2[, or the LEP-LHC ep option. At THERA the interesting x area 
for saturation effects would occur at > 10 GeV^. Another possible way forward is to 



look simultaneously at several processes, for example DIS and inclusive diffraction [RH] or 
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Figure 35: The logarithmic derivatives as a function of compared to the 
predictions of the latest MR ST fits (dotted line) as well as the modification including 
the LO BFKL kernel by Thome (solid line). 

DIS and elastic J/ip production ||129|| . 

Interestingly, the G-B&W model for DIS charm production does predict a turn-over in the 
logarithmic derivative at higher values of Q^. Figure ^ shows a preliminary plot of the 
ZEUS data on DIS charm as a total virtual-photon cross section at constant W vs logQ^. 
There is a clear flattening of the derivative of the cross section in the region somewhat 
less than 10 GeV^, which is well reproduced by the G-B&W prediction. However, this 
effect is not related to the saturation of parton densities at low x, but rather to the charm 
mass and the resultant size of the dipole. In the discussion above on the G-B&W model, 
quark-mass effects were ignored, although they play an important role, particularly at the 
lowest values of Q^, and in ensuring that the cross section matches to the photoproduction 
data. The inclusion of fermion mass effects for instance changes the form of equation |5^ 
to become: 

In 



m}Rl 



where m/ is the fermion mass and rnj » Q^. The effect of the charm quark mass is to 
insert another length scale into the problem ||133| ], since the size of the charm-anticharm 
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Figure 36: The F^^^ data expressed as a total virtual photon-proton cross section 
in bins of constant W as a function of . The lines show the predictions of the 
G-B&W model 



dipoles cannot grow beyond the cut-off imposed by 1/mc, as indicated by equation |9| It 
is this length scale which causes the turn-over in the charm data. Although not related 
to saturation effects, the agreement of the ZEUS data with the G-B&W prediction is a 
beautiful confirmation of the general physics behind the dipole model. 



5 Summary and outlook 



It is difficult to overestimate the effect which the advent of HERA has had on the study 
of low-x physics. It has moved from its infancy at least to and perhaps even beyond the 
kindergarten. One key to this development is of course the vast kinematic range opened 
up by HERA, but the other is the careful experimentation of the ZEUS and HI exper- 
iments as well as the enormous efforts of many theoreticians. The precision of the data 
is now driving many theoretical investigations. Of course, these studies are important 
not only from the point of view of understanding the subtleties of QCD. They are of first 
importance in understand the data from future colliders, in particular LHC Knowledge of 
the data at the kinematic limits of HERA governs the understanding of the backgrounds 
for much of the discovery physics of LHC. The question of deviations from GLAP evolu- 
tion, while fascinating experimentally and theoretically, is also of crucial importance to 
predictions of various SM and exotic processes at LHC. In the next few years, the centre of 
attention at HERA will switch somewhat from low-x to high-x physics, as the HERA up- 
grade allows ZEUS and HI to fulfill their potential as precision probes of the electroweak 
sector. Nevertheless, further precision at HERA and the Tevatron is both possible and 
desirable in the inclusive processes, and the greatly increased luminosity of the upgrade 
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will permit investigations of exclusive processes and difficult areas of phase space which 
will permit the fascinating glimpses of possible deviations from GLAP evolution at low x 
to be investigated further. This is truly an exciting time in low-x physics. 
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